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^> ■ Abstract 

' show that it is possible to derive nontrivial renormahzation group flow from iterative coarse 

ps) ' graining of a closed-time-path action. This renormalization group is different from the usual 

in quantum field theory textbooks, in that it describes nontrivial noise and dissipation. We 
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We investigate the renormalization group approach to nonequilibrium field theory. We 



\ work out a specific example where the variation of the closed-time-path action leads to the 

^ \ so-called Kardar-Parisi-Zhang equation, and show that the renormalization group obtained 

. by coarse graining this action, agrees with the dynamical renormalization group derived by 

O ■ directly coarse graining the equations of motion. 
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1 Introduction 

The goal of this paper is to investigate the renormalization group (RG) approach to nonequilibrium 
field theory. We derive the renormalization group from iterative coarse graining of the Schwinger- 
Keldysh or closed-time-path (CTP) action 0. We work out a specific example where variation 
of the CTP action leads to the so-called Kardar-Parisi-Zhang (KPZ) equation p, ^ ^. We show 
that the renormalization group obtained from the coarse grained action (CGA), agrees with the 
dynamical RG derived by directly coarse graining the equations of motion p, 

The RG [|, 1^, 1^ is a powerful method to analyze complex physical systems. Given a description 
of the system at some scale, a new description at a lower level of resolution is derived by coarse 
graining the former. By analyzing how the picture of the system changes (or fails to change) with 
resolution, important physical information is derived. 

A field theory is most often not considered a fundamental description of a physical system. Its 
field variables are considered as the relevant degrees of freedom at some degree of resolution. This 
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description is not complete, leaving out some uncontrolled sector whose interaction with the field 
variables is characterized as noise and dissipation |I0|. We would like to associate to each level 
of description a corresponding action, so that the changes in this action as we change the resolution 
of our description allow us to define the dynamical RG for the theory. 

A simple way of implementing this idea is by looking at the CTP generating functional, whose 
Legendre transform yields the CTP effective action (EA). The generating functional admits a 
representation as a path integral over fiuctuations in the field variables of the exponential of an 
action functional. By performing a partial integration over some fiuctuations, we obtain a new 
integrand which may be used to define the CGA |[Tl| . The change in the CGA as more fiuctuations 
are integrated away defines the dynamical renormalization group. 

In equilibrium, there is an efficient way to code a description of the system through some 
adequate thermodynamic potential (the free energy for a system in canonical equilibrium, etc.). 
In field theory, the proper thermodynamic potential under canonical equilibrium conditions is the 
Euclidean action, where the time variable is identified with periodicity (3 = 1/T, and T is the 
temperature. The Euclidean CGA is defined from a partial integration over the field variables; the 



variation of the Euclidean CGA with scale is given by the Wegner- Houghton equation |T^, and 
gives rise to the so-called exact renormalization group |13[. 

In dynamical situations, such devices are not forthcoming, and so the dynamical RG is usually 
formulated at the level of the equations of motion 0. Since thermodynamic potentials are most 
often simpler than equations of motion, the equilibrium RG has been much better developed than 
the dynamical RG. 

In dynamical situations, the Lorentzian EA, which may be used, for example, to derive S- 
matrix elements of the field operators, cannot be used to derive a physically sound evolution for 



the background fields . A simple solution lies in adopting the so-called Schwinger-Keldysh 
techniques [|15|, |16[- In this paper, we show that essentially the same ideas can be used to define 
a convenient CTP action for stochastic field theories [0. For the heat diffusion equation near 



equilibrium this was done in [18 



The basic element of the Schwinger-Keldysh or CTP method is the doubling of degrees of free- 
dom. For each field variable in the original theory, a new mirror variable is introduced; accordingly, 
the number of external sources in the generating functional is also doubled, and the EA is defined 
as a Legendre transform with respect to all variables independently. The dynamics for the back- 
ground (also called classical or mean) fields is obtained by taking the variation of the EA, and 
then (but only then) imposing some constraint on the mirror variables, in order to eliminate the 
excess degrees of freedom. The formalism is built in such a way as to make sure that the resulting 
dynamics is causal and respects the reality of the background fields. 

We wish to point out that there are other implementations of the doubling of degrees of freedom 
idea. The best known in this context is possibly the so-called Martin-Rose-Siggia formalism [|l^ 
for stochastic differential equations (SDE), which is closely related to the CTP approach pO |. 

Since the CTP EA may be used to derive real and causal equations of motion for the expectation 
values of field operators (and, in an extension of the formalism, also for their correlations |]10|, ^) it 
is natural to define the RG for nonequilibrium field theory from the iterated coarse graining of the 
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CTP action. This approach to the RG has been put forward in references and [0 . These authors 
show that, under the adiabatic approximation, the RG defined from the CTP action reduces to the 
usual (equihbrium) one. 

However, nonequihbrium field theories also manifest a regime (called strongly dissipative by 
Berera et al. ||2^) with very different properties from equilibrium fields. At the level of the CTP 
EA ||2^, 2|], this regime is characterized by the EA becoming complex, and also by the entanglement 
of the original field variables and their doubles, in such a way that the CTP EA no longer may be 
written as the difference of two independent action functionals. 

These non separable terms are associated to dissipation (when they are real) and noise (when 
they are imaginary). The joint presence of noise and dissipation, which is due to the unitarity 
of the underlying theory, is the dynamical foundation of the fluctuation-dissipation relation near 
equilibrium. 

On closer examination, it is not surprising that studies of the nonequihbrium RG in field theory 
so far have found no evidence of this strongly dissipative regime. The case is analogous to, for 
example, the situation in thermal field theory in 4 space-time dimensions. An approach to the 
renormalization group based on the ultraviolet behavior of correlation functions (which is insensitive 
to temperature ||2^) will fail to disclose the nontrivial fixed point at very high temperatures, when 
the theory becomes effectively three-dimensional. In the same way, the RG derived under the 
adiabatic approximation is insensitive to noise and dissipation, because these are non adiabatic 
effects H. 



In the thermal case, what is needed is an "environmentally friendly" approach to the RG |26 



where temperature dependent correlations and coupling constants are used throughout. In the 
nonequihbrium case, the starting point must be a noisy and dissipative CTP action, including new 
parameters associated to the non adiabatic terms. 

The stumbling block in the completion of this program is the lack of an efficient parametrization 
of the non adiabatic CTP action. Of course it is possible to start from the adiabatic action (that is, 
the difference of two Lorentzian actions for the field and mirror variables, respectively) and derive 
the non adiabatic action after coarse graining some of the quantum fluctuations. However, the fact 
that this is necessarily done in some kind of perturbative scheme (that assumes that the resulting 
corrections to the action are small) to a large extent defeats the purpose of the whole exercise. 
On the other hand, the strongly dissipative regime of nonequihbrium field theory truly exists, to 
such extent that most practical applications of nonequihbrium field theory are actually based in 
stochastic classical field theory, kinetic theory, and even hydrodynamics, all of them limiting cases 
of the strongly dissipative regime. An interesting example is the slow rolling assumption made in 
inflationary cosmology (see, for example, [^). There, the second time derivative of the inflationary 
field, which obeys a Klein-Gordon equation, is discarded when compared with the dissipative term. 

The goal of this paper is to put forward the essential elements of the CTP approach to stochastic 
dynamics, and the derivation of the dynamical RG therefrom. Rather than an abstract presentation, 
we have chosen to work on a specific example. We have chosen a parametrization of the CTP action 
for an scalar field theory whose variation leads to the noisy KPZ equation in 3+1 dimensions P, ^, |^. 
We have chosen this example because it is relatively simple, while its manifold applications warrant 
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its physical cogency |28|, 0. Related with the KPZ equation is the Burger's equation 
pT| which, among its multiple applications, has been useful in describing problems of structure 
formation in cosmology [ p^ . 

The paper is organized as follows. In next section, we introduce the basic notions regarding the 
CTP formalism; then we proceed to defining the CTP action and the CTP generating functional 
for SDE, and show the connection with the usual (single-path) functional formulation for SDE. We 
then apply the CTP approach to the KPZ equation, calculating its associated CTP EA. In Section 

3 we introduce the CG generating functionals and the corresponding actions or CCA. In section 

4 we study the way the CCA runs with changing coarse graining, and compare the resulting RG 
with the one derived by other means. We show that the resulting RG displays nontrivial running 
for the noise and dissipation terms in the action. We conclude with some brief final remarks in 
Section 6. 

In Appendix A it is shown explicitly that the field equations derived from the CTP EA for the 
KPZ equation, reproduce the right dynamics for the classical (i.e. mean or background) field. In 
Appendix B we compute the CGA to second order in the non linearity. In Appendix C we show 
that the effective theory for the modes that survive the coarse graining of the KPZ CTP generating 
functional, is equivalent to that obtained by coarse graining the equations of motion of the field. 



2 Closed Time Path and Stochastic Differential Equations 

After a brief review of CTP formalism, we shall proceed to define the CTP action and the CTP 
generating functional for a class of SDE. Next, we shall apply our method to the specific example 
of the KPZ equation, computing the CTP EA and deriving from it the equation of motion for the 
classical field. 



2.1 Closed Time Path Field Theory 

In the usual In-Out formulation of the Quantum Field Theory, the basic object is the vacuum 
persistence amplitude Z, which encodes all the dynamical information of the theory |Q. Suppose 
we are dealing with a scalar field theory. Then we define 

Z[J] = (Out|In)^ = J V^x) e^5[*l+^^*, (1) 

Here S is the action of the field and J is an external current coupled linearly to the field. (If not 
indicate explicitly, the integrals are over the entire space-time.) This functional generates matrix 
elements of operators between In and Out states, rather than proper expectation values referred to 
a single state. Hence, this formulation is useful when one asks questions about scattering problems 
or rate transitions, for instance. But if we want to deal, in this formalism, with the time evolution 
of true expectation values, we must be able to relate two different complete set of states, e.g. via 
Bogolubov coefficients, in order to relate the In and Out bases. Instead, we can use the functional 



integral method developed by Schwinger and Keldysh, known as CTP formalism [0. In the In-Out 
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formulation, when working in the Heisenberg picture, the vacuum persistence amphtude (0) is also 
given by 



Z[J] = { Out 



Texp 



i j d^x J{x)^Hix) 



In 



(2) 



where T denotes temporal order, and is the field operator in the Heisenberg picture. By 
contrast, in the OTP formulation we define a more symmetric object, namely 



In I In ) 



(3) 



In 



Texp 



dt / d^x J {x)^h{x) 

D J 







Texp 


i f dt 1 




J — OO J 



dt J d^x J+(x)$iy(x) 



In 



That is, one compares the final states that result from the evolution of the In state under the 
infiuence of two external currents, and J~ . Here T means antitemporal order, and t* is some 
late time, which in practice it is chosen to be +oo. It is easily seen that the derivatives of Z[J~^ , J~] 
evaluated at J+ = J~ generate true expectation values of product of fields. In terms of path 
integrals, Z[J^ , J^] has the following representation 



Z[J+, J-] = J V^+{x) V<!>-{x) exp i (^5[<l>+] - S[(!>-] + J J+$+ - J 



(4) 



#+(t*)=$-(t*) 



The quantity Sf^^] — •S'l*!)^] is the CTP action or Sctp- In (§) we integrate over histories and 
$~ that join at time t*. As in the In-Out formalism, the classical equations of motion are obtained 
from the variation of the action with respect to the fields. 
We can define a generating functional 



and classical fields 



-■i\ogZ[J+,J- 



SW[J+, j- 



5J±(x) 

Next, we define the CTP EA as the Legendre transform of W, that is 

n^ti, ^ci] = J-] - 1 + / j-<i> 



d ' 



(5) 



(6) 



(7) 



where it is understood that the currents have been expressed as functions of the classical fields, via 
relations (0). The equations of motion for the classical fields can be written as follows 



tJ^{x). 



The common value of $^ and <l>^, when = J~ , is real, because it is a true expectation value. It 
can be shown from the definitions above that it obeys a causal equation of motion (see for example 
||16||). On the contrary, in the In-Out formalism the classical field is not necessary real, nor the 
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equation of motion that it follows is causal, as it is not a proper expectation value, unless the In 
and Out states coincide. 



It will be convenient to rephrase the CTP formalism in terms of = — $ and (p = <I>^ + $ 
The CTP condition will be given by (pit*) = 0, and the classical equations of motion by 

6ScTp[4>,^] 



6ScTp[4>,^] 



-Jix), (9) 
-J(x), (10) 



where 



Moreover, we shall have 



J^'-^ (111 



SW[J,j] SW[J,j] 
We obtain more symmetrical equations of motion 

l^ch'^cl] = -j, (13) 
T --[(l)cl,^d] = -J- (14) 

When J = 0, the first of these equations gives the physical equation of motion, and the second one 
is trivially satisfied. 

2.2 CTP approach to Stochastic Differential Equations 

The causal and real evolution obtained from the CTP formalism, suggests that a CTP action that 
reproduces the Langevin equation for a stochastic theory described by a field ip, could as well be 
used to compute the correlations of the field and to derive the equation of motion for the classical 
(i.e. mean) field by employing the corresponding CTP EA. In its turn, the defined CTP generating 
functional can be used to implement the RG in the same fashion as it is implemented at the level 
of a thermodynamic partition function for systems in equilibrium. In this way, one would find an 
alternative route for the usual dynamical RG |^. 

Let us consider a stochastic differential equation of the general form, 

r[p]ix)^^-K[p]ix)=vix), (15) 

where K[(p] represents a differential operator not including time derivatives, and where ?7 is a zero 
mean stochastic function with Gaussian probability distribution. It can be seen that this equation 
is obtained from the following CTP action, as explained below, 

ScTp[(l>^v] ~ J dx (j){x)T[(p]{x) + c^- J dx dx' (f){x)N{x, x')(j){x'). (16) 
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Here N{x,x') is the two point correlation function of the noise rj, and c is a dimensional constant 
that makes the action dimensionless. The dimensions of c will depend on the physical interpretation 
we give to (p, e.g. the potential in fluid mechanics, the height function in surface growth. The CTP 
generating functional for this action is 

ZcTp[J,3] = J VcpVip exp iyiScTp[(t),v\ + i J dx [3{x)(f{x) . (17) 

(The CTP condition, 0(t — > oo) = 0, is understood.) To arrive to Eg. ([151) we observe that the term 
exp| — — y dxdx' 0(x)A^(x, x')0(x')| in (p!7D, can be written as the functional Fourier transform 
of an auxiliary functional. We find, up to a constant factor. 



where 



exp \ — — / dxdx' (j){x)N{x,x')(j){x') 



1 



Vrj P[ri] e~*'^/ vix)'l>i^) _ 



P[r]] = exp -j — - / dxdx' rj{x)N ^{x,x')ri{x'] 



(19) 



is the probability distribution of the noise t], and ^ means the inverse matrix of A^. Hence, we 
can write 



ZcTp[3,j] = J VrjVcjyVip P[ri\ exp jiScTpf^, V^, + « j dx [i{x)Lp{x) + j{x)(l){x)] 



where 



CTP\ 



Lp,r]] = c (pix) T[lp]{x) - r]{x) (pix) 



dx. 



(20) 



(21) 



The variation of this action leads to Eq.(|T3]). This method allows us to relate the imaginary part of 
the CTP action, quadratic in the field 0, to stochastic sources [Q, ^ ^, ^ The constant 
c can be absorbed by redefining as c~^0; therefore (p and ip will not have, in general, the same 
dimensions. 

Let us show that the CTP generating functional (|T^) is, up to a Jacobian factor, the generating 
functional usually defined in the theory of SDK from a probabilistic approach fl^, ^ 
namely 



Z[J] = J Vr] P[rj\ exp j dx J{x) (fs{x] rj] 



(22) 



Here ips{x]ri] is the solution -assumed unique- of Eq.(|l5]) for a particular realization of the noise 
7], whose probability distribution is P (19). The derivatives of Z with respect to the external 
current J give the correlation functions of the field, where the average process is referred to the 
noise probability distribution. This formulation of the stochastic problem is equivalent to the 
Martin-Siggia-Rose formalism [|l9l (see ||4^). We now demonstrate that is also equivalent to a CTP 



formulation based on Eg . (|T7|) . Inserting the following identity in (|2^) 

Dip 6\(p{x) - (ps{x;r]] 



(23) 
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we obtain 



Z[J] = J VrjVip P[r]] 5 (p{x) — (ps{x; rj] exp |i J dx J{x)(p{x) 
Changing variables in the argument of the delta functional yields 

Z[J] = J VrfVip P\rj\ 5 T\lp\{x) — r]{x) J exp \i J dx J{x)ip{x 

where 

d SK[ip] \ 



(24) 



(25) 



:r = det|^}=det 



dt 



6ip j 



is the Jacobian associated with the change of variables in the delta functional. If the operator K 
does not contain any time derivative, the Jacobian, up to a field independent factor, is given by 



13,101 



The next step is to expand the delta functional in Fourier components, that is 
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T[^](x) -r]{x] 



and then replace this expression in Eq. 
Eq.(ll), yielding 



(26) 



(27) 



The integral over the noise is done explicitly using 



Vcf) exp <i dx (f){x) T[(p]{x) — t]{x) 



Z[J]= J V(j)V(p J exp l^iScTp[4>,<^] +i J J{x)ip{x) 



(2^ 



where the action Sctp is given by (|I6]) after absorbing the constant c into 0. 

We see that, when the Jacobian J' is field independent, the expression ( p8[ ) can be identified 
with the CTP generating functional we defined before motivated by more heuristic considerations. 
This happens for a broad class of SDE |3^, including the KPZ and Navier-Stokes equation 

H- 

It can be seen that eqs.(||) and (^) are equivalent to the equations proposed for the physical 
and the non-physical field operators, respectively, in the work of Martin, Siggia and Rose To 
show this, we adopt the notation of that paper, thus 

K[ip]{xi) = / dX2 U2{xi, X2)(p{x2) + / dX2dX3U3{xi,X2,X3)(p{x2)(p{x3), 



(29) 

and Ui{xi) = r]{xi). Hence, remembering that the noise term in Eq.(^) is recovered after Fourier 
transform the quadratic term in (j), the classical and physical (i.e. J = 0) equations of motion 
derived from the CTP action -eqs.(|D and (^), will be 



6S, 



CTP 



S(j){xi) 
5 Sctp 
5ip{xi) 



9^(xi) 

dt 
d(t){xi) 

dt 



dX2U2{Xi,X2)(p{x2) - J dX2dx3U3{Xi,X2,X3)(p{x2)(p{x3) = f/i(xi),(30) 

- / dx2 (l)ix2)U2{x2, xi) -2 / dx2 dx3 U3{x2, X3, Xi) 0(x2)v?(a;3) = 0. (31) 



8 



These are the same as eqs.(2.1) and (3.1b) in ref. 



To conclude this section we compare the CTP approach with the single-path approach to SDE 
of (see also ^ 0). The difference arises in Eq.(^). If the integral over the noise is 



performed explicitly with the aid of the delta functional, we get 

Z[J] = JVip P[T[ip]\ J expji J dx J(x)(/?(x)| . (32) 

Using the definition ([T9|), and leaving apart the Jacobian J', we obtain the following single-path 
action 

Ssp = -\j dxdx' T[ip]{x) N-\x,x') T[^]{x'). (33) 

While it is a valid representation of the generating functional, this action can not be used to 
generate the dynamics of the classical fields. To see this, suppose that the operator K in (plSj) is 
linear in ip, and moreover that the Green function G, associate with that equation, is causal. The 
assumption regarding the linearity of K implies that S'5p[(y9] is quadratic, and that the Jacobian 
J' is field independent, so it can be ignored. Because T is a hnear operator, the classical field 
associated with Ssp will obey the classical equation of motion obtained from the variation of Ssp- 
So, when an external current J is coupled to the field we obtain 

M^) = G(.,x0iV(..,..)Gt(...X3)J(.3). (34) 



Because of causality the Green function verifies G{x,Xi) oc 9(t — ti) and G\x2,Xs) oc Oit^ — t2), 
and hence we can not affirm that if obeys a causal evolution. Therefore the physical meaning of ipci 
is limited, as for the classical fields in the In-Out formalism of quantum field theory. This problem 
is related to the fact that operator T appears twice in (^). As noted in [0, the set of solutions 



associated with the variation of Ssp in (|33|), includes not only the solutions of (p!5|), which are 
causal, but a set of spurious solutions. 



3 CTP Approach to the KPZ Equation 

In this section we apply CTP methods to the KPZ equation 

^-.VV-^(Vv.)^ = r^, (35) 

where t] is the noise term, assumed to have some particular Gaussian statistics. The KPZ is a 
well known equation that belongs to a general class of stochastic non-linear differential equations 
of diffusive type, for which our method can be extended straightforwardly. 

There is a large amount of literature regarding the KPZ equation (see, for example |^ and 
references therein). We mention here only a few points concerning it. In the context of fluid 
dynamics, the KPZ equation is derived for the case of free-vorticity, null-pressure fluid, — Vy? being 
the velocity fleld. When used to describe some phenomena related to surface growth 0, |[, the 
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KPZ equation is also derived as one of the simplest non-linear extension of the Edwards- Wilkinson 
equation, (p measuring surface height. In addition, the KPZ equation is closely related with flame- 
front propagation |2^, dissipative transport and polymer physics |^, to quote some. When 



derived from Navier-Stokes equation, it is seen that the non-linear coupling must be A = 1, which 
is not necessarily the case in treating, for example, surface growth. The noiseless KPZ equation 
is Galilean invariant, a property which, in the case of a fluid, is inherited from the Navier-Stokes 
equation, and that in the context of surface growth is related to the rotation symmetry of the 
coordinate system. If the noise is white, and translation invariant, this symmetry is preserved by 
the noisy KPZ equation as well. This fact implies a non-perturbative result concerning the running 
of the coupling A when the RG is implemented 0, thus reducing the number of independent 



scaling exponents. As in the paper of Forster, Nelson and Stephen for the case of Navier-Stokes 
equation |Q (see also one can derive the RG equations by coarse graining the equation of 



motion of the field (p. We want to apply ideas concerning the CTP formulation of quantum field 
theory in order to implement this RG transformation at the level of a CTP generating functional. 

We begin defining the CTP action and generating functional for the KPZ equation. The free 
case is examined in order to implement the perturbative calculation of the EA for the interacting 
case. 

3.1 CTP Action and Generating Functional for KPZ Equation 



As demonstrated in the Jacobian J associated to the KPZ equation is field independent, and 
hence we can define a CTP action for the KPZ equation which describes the stochastic dynamics 
of the field ip. In 3+1 dimensions, absorbing c into Eq.(|16D yields 

^CTp[0,<^] = / d^x \^4^{x)C^{x) - ^</)(a;)(V<^)2(x)| + d''xd^x'(j){x)N{x,x')(j){x'), (36) 

where 

£ = |-.V^ (37) 

N{x, x') is the two point correlation function of the noise rj, assumed Gaussian and having zero 
mean value. With these definitions and per the early discussion we see that the KPZ equation is 
attained. 

3.2 The free case 

If the non-linearity is absent, i.e. if A = 0, we are dealing with the free case, and the corresponding 
free action is 

So[(f), if] = J d^x (f){x)C^{x) + d^xd^x' (t){x)N{x, x')0(x'). (38) 

When linear couplings of the fields with external currents are included, the variation of the free 
action with respect to the fields, gives the classical equations of motion, namely 

5Sq 



50(x) 



C^{x) + / d^xN{x,x)(t){x) = -j{x), (39) 
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where 

^* = Ql + (41) 



We can write the solutions to eqs.(|39|) and ( ^0|) in terms of the fundamental solutions G and G*, 
which satisfy 

L^G{x,x') =LlG*{x,x') =5^(x-x'). (42) 
Explicitly (in 1+3 dimensions) we have 



Hence, 

0(a;) = f d^xiG*{x,xi)3{xi), (45) 



(There is the freedom on add some arbitrary solutions of the homogeneous equations. However, 
the only such solution which is bounded for all times is identically zero.) 
If J = then = 0, and (p is given by 

= (J.|n!|)3/^ ^(t-Mjfa). (47) 

This entails a causal evolution. 

We define the generating functional Zo[J,j] for the free fields 

Zo[3,j] = [ V(j){x)Vip{x) e^5o[<^'^]+^/(j'<^+-i^) . (48) 



The mean fields are obtained by differentiating —i log Zq with respect to the currents. Having in 
mind that for the free case the mean fields satisfy the same equations (|39|) and (0), we find (up 
to a normalization factor) 

Zo[J,i] = exp ^ y rf^xirf^X2< - J(a;i)G(xi,X2)j(x2) +i(xi)G*(xi,X2)J(x2) 



iJ{xi] 



d X^d X4G(Xi, X3)A^(X3, X4)G*(X4, X2) 



3{X2)\. (49) 
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The two point correlation functions are given by the second derivatives of Zq, so we have 

{(p{x)(p{x')) = — d'^xi / d'^X2G{x, Xi)N{xi, X2)G*{x2, x'), (50) 



{(f){x)ip{x')) = —iG*{x,x'). 
The corresponding functions in the momentum space are 



{<l>ipMp')) 



[We indicate the Fourier transformed fields with the same name as the original fields, and use the fol- 
lowing convention in d+1 dimensions: f{p) = (27r)^*^'^+^)/^ J d'^x dx^ exp |— i (j)^x^ — p ■ xjj f{x) , 
where x^ = t.] 

3.3 The interacting case 

Let us go back to the definition of the generating functional for interacting fields (dropping the 
CTP subscripts) 

Z[J^j] = J T)(f)V(p e*'^[''^''^]+* / '^'^^ [K'-^)H^)+^{x)'p{x)] ^ ^^52) 

The EA is given by 

rife. = log Z - / A + ^.(x)J(.)| . (53) 

and admits the CTP representation 

^n<P.i,Vci] ^ f i)0(a;)r'(^(x)e*^['^^'+'^' '^='+^1 , (54) 
Jipi 

where IPI indicates that only diagrams one particle irreducible must be included in the diagram- 
matic evaluation of the functional integrals. 
For the KPZ CTP action (BBI) we have 



A /■ r 

S[(j)ci + 4>, ^ci + ¥^]= 5'[0d, V^d] + So[(f), ip] -1^ j 'P {^^f + 'Pel (y^f + 20V ■ Vipd 

+ linear terms in and if . (55) 

Taking the logarithm of the Eg. ([5^) and expanding to (9(A)^, it results 

A f 

r[(f),i, ifci] = S[(Pci, fci] --Jd^x{[(f) {V<ff + (Pel + 2Vif ■ V<fci)^) 

' d^xd^x'{{<j) {Vipf + (j)ei {Vipf + 2<f)Vip ■ V^ei)^ 

X (0 (Vipf + (Pel (V^)' + 20Vy^ ■ ^^ci)^,)^^^^^^^^ , (56) 
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+ 



where the averaging operation (...) is defined as 

Note that the term ((^Vv? ■ Vv?d)x(0Vv2 ■ VvJd)^')' which could give a nontrivial equation of 
motion for (pcu vanishes, because it is proportional (up to spatial derivatives) to 

(0(x)v?(x')) (0(x')¥5(x)) oc e{t' - t)e(t - 1'). 

Since 5*0 is quadratic, expectation values may be written as products of the two-point correlation 



functions given in ( ^OD and (|5T1) . Hence, after Fourier transformation of (pq), the result is 

r[0d, Vd] = s[(j)ci, (fid] + ^s[(f)d, ifid] , (5^ 

where 

A 



AS[(j)d,^d] = ^ d^pid^p2Aii{Pi -P2,P2) (pdi-Pi) 



I / /"^4 >4 -A^j(p2, -P2) {Pl)i{Pl+P2) 



A? 

+ Up^] [-{pi + P2f + y{pi + P2)^] 
,4 >4 >4 ^ii(P2,P3) (P2 -Pl)i(-Pl +P2 +P3)i 



+ Ai d pid p2d p3^ — — — — —-^ -(l)d{-Pi>^d{Pi -P2- P3> 

J \[pi - P2r + i^[pi - P2Y\ 

+ j d^Pid^p2d'^P3d'^P4Aij{p2,P3) ^ij {Pl - P2,P3 - P4:)4>d{-Pi)4>d{-P3) |- (59) 

The sum over repeated indices is understood, and 

^-(^'^'^ = [.p0 + ,^i]|^p'^i + .p^2] M • (60) 

The equations of motion for the classical fields result from the first variations of the EA. Those 
with proper physical meaning are obtained when 0, J = 0. Thus, to (9(A^), 

-[(pel = 0,(Pcl\ = 



Hdi-p) 

X p X f 

[ip^ + iyp^](Pd + J d^pi Pl ■ {p - Pl) y^diPi) Vd{p - Pl) + J d^piAii{p-pi,pi) 

/■,4 .4 ^ij{PuP2){Pl-p)t{Pl+P2-p)j , \ ■( \ fai\ 

' d pid P2 TT^- — —. — — ^d{p-Pi-P2) = -j{p)- (61) 



IGtt^ 7 [{p — pi)'^ + ^{p — pi^"^ 

The equation 



5T 

' d = 0, ^d] = (62) 



5Vci{.-p) 
is automatically satisfied. 

In Appendix A we show explicitly that Eq. (|6lD is also obtained by averaging the KPZ equation 
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4 Coarse Grained CTP action for the KPZ equation 



In order to implement the renormalization group (RG) transform, we analyze the influence that the 
modes of higher wave number exert on lower ones, by computing the Coarse Grained Action (CGA) 
|T], 1^, [IT], When we are only concerned with the lower wave number sector of the theory, we can 
carry out explicitly, in the generating functional Z, the integration over the higher wave number 
modes, and the result of this partial integration will be a functional of the lower wave number 
modes only. This functional is indeed a generating functional for the lower modes, in which the 
influence of the higher modes is incorporated as modifications of the original action. 

This procedure may be seen as a straightforward application of the Feynman- Vernon influ- 
ence functional techniques to this problem [^, where the low wave number sector is regarded as 
"system" and the short wave modes as "bath". 

To the best of our knowledge this approach has not been systematically discussed in the lit- 
erature on SDE. There exist works where the RG is also derived from functional formulations of 



the stochastic theory (see, for example, |^ for critical dynamics of Helium, antiferromagnetics 
and Liquid-gas systems; or |Q for the KPZ equation). The crucial difference between these 
works and the present paper, is that we coarse grain the generating functional explicitly, imposing 
an ultraviolet cutoff, while in the cited papers the RG is obtained from the study of the singular 
ultraviolet contributions to the many-points response functions. 



We start with our early definition ( [5^ ) of the generating functional for the interacting fields 

Z[J,j] = J V(j)V(p e*^['^''^]+*/'^*'^ \j{x)^{x)+3{xMx)]_ (^g3^ 
Now we split up the field and currents, according a scale that we shall choose below, 

= 4>> + 4><, 

3 = J>+J<, 
J = J> + J<. 

Here, ipy contains the modes of higher wave number, contains the lower ones, and analogously 
for the other quantities. The division will be specified by a cutoff 

^<(a:) = r 4^,e^^'^'^'-^'-'^^{k',k) (64) 
J\k\<As (27r)^ 

^>{x) = [ ^^,e^^''^'--'^^{k\k) (65) 

and so on. Here, i~p{k) is the Fourier transform of the field '■p{x), and A can be identified with a 
natural cutoff of the theory. 

In any case, the correlations are obtained from the variation of Z with respect to the currents, 
and after that, by setting the currents equal to zero. As stated early, we just want to compute 
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correlation functions involving the lower wave number modes. That can be accomplished merely 
by the variation of Z with respect to j< and J<. Therefore, it will be enough if we set, from the 
beginning, j>, J> = 0. The CGA is achieved by performing explicitly the functional integrations 
over and 0>. We rewrite the action 5'[0, = 5'[0> + 0<, + yj^] in the following manageable 
way, which will be useful to compute the CGA perturbatively, 

S\<\>> + 0<, V5> + V5<] = 5'[0<, V9<] + S'o[0>, V5>] + 5'/[0>, 0<, v?>, V9<]. (66) 

Here, corresponds to the free action of the original theory. Hence, it results 

^[J<,J<] = 

Vip^V(f)^ g*'S'[9i<, ■/'<]+« / d'^x\j<<t><+J<ip<c]{x) ^ T>(p~^V(j)-^ ^iSo[<t>>,<fi>]+iSi[<t><,<l)>,ip>,tp<] 

= j V(^^V(j)^ ^iS[4><:,<^<,]+i J ci'*x[j<</><+J<v?<](x) ^ giA5[</>< , i/3<] ^ ^g^) 

The CGA is defined as 

5'cg[0<, = S[(f)<, + AS'[0<, v?<]. (68) 

In the present paper we are concerned with the KPZ equation and with its associated CTP action 
(|36|) , which in p-space is given by 

S[(f), Lf] = J d^p (l){-p) (ip^ + z/p 2) Lp{p) 
A /■ 

+ 2 J d'^P^d^P'^d'^Ps (27r)"^5(pi + P2 + Ps) P2 ■ Ps 0(Pi) v{P2) viPs) 
+ d*p,d% 0(pi)Ar(-pi,-p2) 0(P2). (69) 



2 

Splitting the fields according to wave number yields 



S[(f), (f] = J d^p (t)>{-p) (ip^ + z/p ^) v3>(p) + j d'^p (t)<{-p) (ip^ + z/p ^) v2<(p) 

+ \j d^Pid'^P2 <P>{Pi)N{-pi, -P2)0>(P2) + ^ / d'^pid'^p2 0<(pi)iV(-pi, -P2)0<(P2) 



+ i j d%d% (f)y{pi)N{-pi,-p2)(l)<{p2) 

J d'^Pld'^P2d% (27r)"^5(pi +P2 +P3) ■ P3{0>lV5>2V5>3 + 0<lV5>2V5>3 + 



A 
2 



+ 2(/)<iV9>2V2<3 + 20>iV?>2V5<3 + 0>2V5<2V5<3 + '/'<lV5<2V5<3}- (70) 

We shall assume that the noise is translation invariant (TI), therefore the term in the third line is 
zero because of orthogonality. Hence, as before, we have 

5o[0>, <^>] = I d^p 0>(-p) (zp° + z/p 2) + ^ 1 rfV^V 0>(pi)iV (-P1, -P2) 0>(P2), (71) 

and using the definition given in Eg. (p6|) , we find 

A f 

Si[(t)>,<P<,^>,^<] = J d'^pid^d'^ps (27r)"^5(pi +P2 + P3) P2 •P3{0>i¥^>2¥'>3 

+ 0<lV5>2V5>3 + 20<lV5>2V5<3 + 20>lV5>2V5<3 + 0>2V5<2V5<3} • (^2) 
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Therefore, from Eg. (]67|) , we obtain 

giA5[<^<,v:)<] ^ J X>(j)^Vip^ gi5o[</.>,</p>] 

X exp 



i\ f 

— J (fpid'^p2d'^Pz {2txY'^5{pi+P2+Pz) P2 ■ Ps {</>>i</?>2<^>3 

When the noise is white, TI and have no spatial correlations, we have 

N{p,p') = 2D 5{p'' +p"') S{p + p'), 
where D is the noise amplitude. For this case (see details in Appendix B): 

AS[<j>^,^^] = -^Jd'p{27r)'S{p)J'Mp)+'-Jd'p 0<(-p) 2SD 0<(p) 
+2D\^ f d^p <p<{-p) ip<ip) X p^5u{p) 



X f 

— Y J d'^pd'^qd'^kdH C{p, q, k,l) (j)<{-p) ip<{q) ip<{k) ip^{l). 
We have defined the tadpole amplitude, 

d^q 2Dq^ M{q) 



the noise amplitude correction 



167r4 [(g0)2 + ^2(^2)2]: 



2An = 9n2),2 / A [g-(p-g)]' M(g,p-g) 

7 4712 [(^0)2 + ^2(^2)2] [(p0_g0)2 + ^2{(p'_^2}2]^ 

the function A, related to the arising of multiplicative noise, 
the viscosity correction 

d^g q-{j) — q) p ■ {p — q) ^{q,p — q) 



p 8v{p) 



167r4 (igO + zyg 2)[(p0 _ ^0)2 + ^,2|(p'_ ^2}2] ' 

and the (/^^-interaction coupling 

q- (k + l) k-l yi{k + l) 



C{p, g, k, I) = (27r)-^ S{-p + q + k + l) 



[i{ko + P) + u{k + If 
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In its turn, M{p, q, . . . ,k) means that the momenta in the set {p,q, . . . , k} are restricted (i.e, must 
be projected) to the momentum shells < |p| < A, < |g| < A, and so on. 

In conclusion, when ( |63[ ) is coarse grained, the generating functional for the remaining modes 
is obtained by modifying the original viscosity u and the noise amplitude D, and by adding some 
new terms: a tadpole term that concerns the homogenous mode only, a multiplicative noise term 
(see Appendix B), and a cubic interaction term. We remark that the noise terms are read directly 
from the imaginary part of the CTP CGA. 

5 Renormalization Group from CTP CGA 

The action S we started with (P^), is actually a coarse grained action. The fields ip and are 
assumed to describe the physical world up to certain degree of resolution, hmited, eventually, by 
a natural cutoff A, as can be the atomic size in a turbulent fluid or the Compton length of heavy 
particles in particle physics. When we integrated the higher wave number modes in the generating 
functional, we obtained a new action, suitable for a physical description with a lower degree of 
resolution A^. 

Suppose that we are interested in the behavior of the theory at momentum scales not superior 
than e~* A, with s real and positive. In principle, this implies that, in the integrations we performed 
in the Section ^, some linear combinations of the momenta must be restricted to the shell e~* A < 
IpI < A. Often we are only concerned with the small-p modes, for which p is near to 0. Hence, we 
must integrate all the modes except those very close to the origin, as close as necessary to obtain 
a leading order result. However, it can be seen that in the case of the KPZ equation, as in others, 
divergences arise in the limit of A^ —>■ 0, indicating that the perturbative approach fails p, ^ 

What we can do instead, is to implement the so called RG formalism 0, 0, ^. The scheme is 

(i) to perform the integration from A to A not at once, but in repeated integrations over 
infinitesimal shells in 3-momentum space. In integrating over one such shell, the cutoff changes 
from A' to e'^A' ^ {1 - 6s)A'. 

(ii) After each shell is integrated, the fields, lengths, times, momenta, etc., must be rescaled to 
bring the theory to its original aspect. In particular the rescaling of momenta must adjust the 
cutoff to its initial value at the beginning of the process, and, in addition, some factors can affect 
the coupling constants. 

When combined and repeated these operations give sensible results. It must be clear that we 
are not simply calculating an integral as the sum of discrete contributions from a partition of the 
domain, because at each step the coupling constant that measures the perturbation is renormalized, 
that is, at each step we are perturbating with respect to a different coupling constant; it is an 
iterative process that gives meaning to the whole integration between A and A — 0. 

In section ^ we have already performed the first step of the RG scheme: the coarse graining 
of the generating functional. The result was that the generating functional for the long modes is 
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obtained from the original one by introducing i) a correction to the noise correlation function, given 
in Eq.(|7^) ii) a correction to the viscosity- like coupling, given in Eq.(|79D, and finally by including a 
set of new terms: the first (tadpole term), the third (multiplicative noise term) and the fifth (cubic 
interaction term) in Eg. ([75|) . Extra terms are a common by-product when one coarse grains a 
generating functional . Actually, the tadpole term can be eliminated by a simple transformation, 
V<{p) ~^ ¥^<{p) + ^27r^A JF(9o5(p); thus only the multiphcative noise (MN) and the cubic interaction 
(CI) terms remain. At this point, if we proceed further and repeat the coarse graining, it can be 
seen that, at (9(A^), no others terms arise. The effect of the MN and that of the CI is just to correct 
the terms already present in the CCA in a way that can be traced systematically. Hence, the first 
time we do the coarse graining is very special, because the effective viscosity is now a momentum 
dependent function, new terms arise that were not included in the original action, and no other 
terms appear when we repeat the coarse graining. 

The natural question is why not to include this momentum dependence and the new terms from 
the very beginning. The momentum dependence of the viscosity can be ignored if one is interested 
in the p — > limit of the theory only. On the other hand, the MN and the CI terms, because of 
the constraint that some momenta must be on the shell, involve modes that despite being <, must 
lie close to the shell. Hence, if it is assumed that the fields < have support near the origin, these 
extra terms vanish. Moreover, under certain assumptions MN and CI terms, as more and more 
shells are integrated and the variables rescaled, tend to vanish, i.e. they are irrelevant terms. We 
shall not include them in our treatment of the RG (see below). 

In what follows we shall work out in an arbitrary number of spatial dimensions d, and, further- 
more, assume that initially A = 1, with the appropriate dimensions. As before, the noise verifies 



Let us study the small-p limit of the corrections introduced by the coarse graining. We start 
with Eg. ([77|) . There, q and p — q must be on the shell between (1 — 6s) A and A, and because we 
shall restrict ourselves to the small-p limit, we must inspect the behavior of the integral when the 
external momenta p get close to 0. To lowest order, the effective noise satisfies ([7^) provided D is 
adjusted by the following amount [Q, |, |] 

5D = ^^6s, (81) 

where Kd = Sd/i^irY, and Sd is the area of a unit sphere in d dimensions. A similar conclusion is 
reached for Eq.([79|) ^, H; in the small-p limit, we find that z/ must be replaced hj u + Su, where 

= -Kd^^Ss. (82) 
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Hence, when attention is paid to the small-p modes, the CGA will be given by (eliminating the 
tadpole, discarding MN and CI terms and dropping the subscripts <) 

Scg[^: v] - I d'^'p H-pW + W + H p'Mp) 

j d'^^^Pld'^^^P2d'^^^P3 S{pi+P2+P3) P2-P3 0(P1)<^(P2)<^(P3) 

d'^+V [2D + 25D] (t){p). (83) 



A 

i 



+ 
+ 

Next, we proceed with the rescaling. We take 6 = 1 + (5s, and define 

Vip\p) = h^^^+'^ip'^p'), 

<l>ip',p) = r"+^0(p'°,p'), 

where = p° and p' — b p. Therefore, we obtain 

Scg[(I), (p] = S[4), (fi] = 

' d'^+^p 4){p) {ip^ + {b'-^ [u + 5u]} (^(p) 
A 



(84) 
(85) 



I' 



+ 



87r2 



J d'^^^pid'^^^P2d'^^^P3 6{pi +P2+ Ps) P2 ■ P3 <j){Pl) ^{P2) ^{Ps) 

J d'+'p 0(-p) {b-'^-'+^[2D + 25D]} ^(p). (86) 

Some remarks are in order. The new variable p is such that \p\ runs up to A, as for the original fields. 
The exponent a + z + d, that rescales the field (p{p) in the momentum representation, matches with 
an exponent equal to a for the rescaling of (p{x). The choice of the exponent given for 0, makes 
the free part of the rescaled CGA form-invariant, and hence, we can iterate the process without 
further modifications. 

In conclusion, after integration and rescaling are performed, the action is characterized by a 
different viscosity 

j) = {1 + 6sy-^ [i^ + 6i^], (87) 

by a new coupling constant 

A = (1 + (5s)"+^-' A, (88) 



and by a new noise correlation coefficient 

L> = (1 + (5s) -2"-'^+^ [D + SD]. 



(89) 



These are general relations, which are valid for every two consecutive instances of the RG procedure. 
Finally, we can arrive to a set differential equations for the running of these quantities, namely 



du 
ds 



V 



Z-2-K, 



X^Dd-2 



dX 

— = X[a + z-2], 

as 



dD 
ds 



D 



z — d — 2a + 



Ad 



x^dk; 



(90) 

(91) 
(92) 
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This are the well know RG equations for the KPZ equation 0, |^. 



In the analysis we made above, we discarded some terms that result from coarse graining the 
initial generating functional for the KPZ equation, the MN and the CI terms. In principle it is not 
difficult to take into account their effect in a systematical manner (for the Navier-Stokes equation 
see |^9[). On the other hand, provided ci > 2 we can see that that both the MN term and the CI 
are irrelevant in the special case when we are near the trivial fixed point. This fixed point is given 
by A = 0, 2 = 2 and a = 1 - d/2. The MN term rescales as h^'^-'^\ and the CI as h^'^~^'^\ Hence, if 
d is greater than 2 both terms tends to zero exponentially when s ^ oo. 



6 Final Remarks 

In this paper we accomplished two things: 

a) with respect to the theory of the nonequilibrium renormalization group, we show that it is 
possible to derive a nontrivial renormalization group fiow from a CTP action. This renormalization 
group is different from the usual in quantum field theory textbooks (see for example |5^) in that 
it describes nontrivial noise and dissipation. This regime has not been observed in earlier studies 
of the renormalization group from the CTP effective action []I| . In these studies, the starting point 
was a noiseless, time reversal invariant theory, which was investigated within perturbation theory. 
But the relevant noise and dissipation effects are essentially nonperturbative 0]. A nontrivial 
nonequilibrium renormalization group can only be found in an "environmentally friendly" approach 
[ p6| where the basic description of the theory already has noise and dissipation built in. 

b) From the point of view of the renormalization group flow in the KPZ equation, we have 
derived the relevant flow equations from an analysis which consistently considered only the long 
wavelength sector of the theory. The usual approach of deriving these equations from the ultraviolet 
behavior of response functions although technically correct, is conceptually contrived. In this 



approach, the renormalization group flow is derived from a regime where the noisy and dissipative 
effective description embodied in the KPZ equation ceases to be valid, and the underlying "unitary" 
theory is recovered. On the other hand, the ultraviolet divergences in this underlying theory ought 



to be the same as in vacuum, therefore leading to the usual "textbook" renormalization group ||50|] . 
For this reason we believe the approach in this paper, where no reference to ultraviolet behavior is 
made, is conceptually simpler, although technically equivalent. 

One thing we did not accomplish is to describe in detail the crossover from the high-energy 
unitary theory to the low-energy noisy and dissipative effective theory. We bypassed this difficult 
problem by choosing as low-energy effective description a theory with a clear physical content. For 
example, if the high-energy theory leads to hydrodynamics in some regime, then it contains the 
Burgers and KPZ equation in the limit in which the pressure is null and the velocity field vorticity 
free. If we consider a the theory of a scalar field, for example, we know this limit exists, because 
at low temperatures the field will behave as a condensate and develop a negative pressure, while 
at high temperatures the theory will be approximately conformally invariant, thus leading to a 
radiation-like equation of state. Thus the pressure will be much lower than the energy density at 
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least in some intermediate range. The KPZ field, of course, is a collective mode when described in 
terms of the fundamental theory. We expect to continue our research on this issue. 

The renormalization group as studied in this paper is a necessary tool to understand the nature 
of collective variables describing the relevant physics in strongly interacting nonequilibrium systems 
such as the Universe during the reheating period and the gluon fireball in the early stages of a high- 
energy heavy ion collision. We continue our research on this rewarding problem. 
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8 Appendix A 

8.1 Average of the Langevin KPZ equation 



As a way of comparison with the results of Section p.3| , in this appendix we shall calculate the 
equation of motion for the mean (i.e. classical) field directly from the noisy KPZ equation (|35|). In 
momentum space it reads 



A 



rfV Pi ■ {p- Pi)'f{pi) v{p - Pi) = vip)- 



(93) 



We write ip for the mean value of (p after average out the noise rj, and define the fluctuating field 
ip according to (f = lp + ip. Therefore, if we average the KPZ equation, it yields 



v>{p) 



A 



And thus 



i){p) + 



A 



(fpi pi-{p- Pi) [p{pi)^{p - Pi) + (V'(Pi)V^(p - Pi))] = 0. (94) 



d^Pi Pi ■ {p- Pi) [2tp{pi)ip{p - pi) (95) 
+ ^(Pi)^(p - Pi) - {^{Pi)^{p - Pi)) ] = V{p)- 



We shall write the solution for ijj as a. power serie in A, 

^(p) = ^(o)(p) + A^(i)(p) + 
from which the following expressions result 

^pW^p) 



(96) 



[ip'^ + up'^] 
1 



(97) 



47r2 [ipO + vp' 



d^Pi Pi ■ (P-Pi) 2V'(°Hpi)^(p -Pi) + ij^''>{pi)ij^''^ip-pi) 



(0)/ 



i^^'\pi)^''>ip-pi) 



(0), 



(9^ 
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Hence, 



Back to Eq.(pl) we find that 



N{p,p') 



[ip^ + up'^] [ip' + up' ^] ' 
An^[tp^ + up^]J "^^'^^ ^ 



[ipi + upi] [ip'^ + up' 2] 



(99) 
(100) 



+ up"^ 



^(P) + ^ I ^^Pi ^1 ■ ~ ^i) ^(Pi) ~ Pi) + ^ / d^Pi^ijiPi,p- Pi)k 



,4 ,4 Aii(P2,P-Pl) (Pl -P2)iPl, 

" Pid P2 



[ipi + up^ 



^{Pi -P2) = 0, 



(101) 



where A is defined in (|60|) . The change pi —pi + p yields to the same expression we found by 
computing the CTP EA, Eq. (|6lD , when j = 0. 

9 Appendix B 

In this appendix we evaluate Eq. (|T3D to order and compare the resulting CG A with the equations 
of motion obtained from coarse graining the KPZ equation. 



We start from Eq. ([73|) . The average of an odd number of > fields is zero, because of the parity 
of the free action (^) when we change V9> by — and, simultaneously, 0> by — 0>. Hence, up to 
C>(A2), we find 



A 



i5o[</>>,v> 



X {0>1 V5>2 V>2. 0>i V5>2 V5>3 

+ 0<1 0<i ^>2 V5>3 V5>2 V5>3 

+ 40<i V?<2 0<1 V5<2 V5>3 V5>3 

+ 4V9<3 V5<3 0>i V5>2 0>1 </'>2 

+ V5<2 V5<3 V<2 ^<3 'Pyi 'Pyi 

+ 40<i V5<2 0>i V3>2 V5>3 V5>3 

+ 2V9<2 V5<3 0>i V5>2 V5>3 0>1 

+ 40<i V9<3 V5>2 V5>3 0>i (/'>2 

+ 4 0<i (y9<2 V?<2 V'<3 V5>3 0>i} 



C?Pl23 P2 ■ P3 {4><lf>2'^>3 + 20>l</'>2</'<3} 
y" dpi23 dqi23 P2 ■ P3 gs ■ q3 



(102) 



connected 
j4^ M„ 



Here, (v5<i) means (p<{pi) i^<iqi))', dpus stands for d pi d p2 d p3 5(pi+P2+P3) and analogously 
for dqi23- In computing the last expression, as the logarithm of (|73|) is taken, we must discard the 
disconnected diagrams associated with each functional integration. Also, we must have in mind 
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Pi 



(2ji)-2 X S(pj + Pj + p^ ) PtP^ 
2 ' ' 



P3 



:up)%{p')}= ---- 



>>(i?) 9>{pX 



Figure 1: Vertex and propagators used to calculate the coarse grained action for the KPZ equation. 

that when some integration variable, say pi, satisfies < A^, then ^>{pi) and </>>(pi) will vanish; 
conversely if > A^. To evaluate (|102|) we can employ the propagators given in ( ^I]) with slightly 
modifications, that is 

(0>(p)0>(p')) = 0, 
{<P>{P)V>{P')) 



(27r)2[ij9'0 + z/p'2]' 

It is understood that if some momentum in the previous equations lies below A^, the corresponding 
propagator will be null. In Fig. 1 we give the convention adopted to represent the propagators 
listed above. These propagators will be used as internal lines in Feynman diagrams. On the other 
hand, when computing these diagrams, for each vertex there will be an integral over the three 
momenta attached to it. After splitting the fields, the propagators quoted in (|50| ) are not valid any 
longer. However, because the split is in wave length and not in frequency, the causal properties of 
the propagators are still the same. 



Consider the terms of order A in ( |102|) , the first of which adds to the action a term that is 
(functionally) linear in 0<, 

(2.)-^ Hp) MP) X [I^JS "^'^(o^' . (1°^) 

where we have assumed that the noise represented by N is not only TI but also white. Diagram- 
matically this term is shown in Fig. 2a. The external lines take trace of the < fields that are 
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(a) 




Figure 2: Order A Feynman diagrams for the coarse grained action. The external fields are indicated by 
double lines. 

attached to a given vertex: a double continuous line represents a field; a double dashed line is 
used for a 0< field. The contribution given in ( |104| ), in turn, can be seen as a field-independent 
term added to the classical KPZ equation (0). The remaining term of 0{X), represented in Fig. 
2b., is zero. This is because the propagator {(p>{pi)f>{p2)) introduces a 6{pi +^2)- In addition 
with the conservation delta, 6{pi +p2 +P3), it implies a S^ps). The product of this S{p2,) with p2 -ps 
force the whole integral to vanish. We now proceed to consider the terms of order in Eq.( |102[ ): 

a) The first one corresponds to three connected, non equivalent diagrams (Fig. 3a.), and gives a 
contribution does not depend neither on nor 0<. However, each of these diagrams is zero, either 
because they entail a delta evaluated in a momentum that lies outside the integration domain, or 
a product of two mutually excluding ^'s. We shall find more of these cancellations below. 

b) The second term of order A^ (Fig. 3b.), consists on a closed loop, and has the same structure 
as the noise term in the original action. Explicitly, this diagram gives the following contribution to 
the CGA action 

^y"rfVc?^?i 0<i0<i d'^P2d'^P3d^q2d'^q3 (27r)"''5i23 5i23 P2 ■ Ps q2 ■ Qs 

N{p2,q2)N{p3,q3) \ 

where ^123 stands for S{pi +P2 + P3) and 6123 for 6{qi + q2+q3)- M is the product of the projectors 
over each 3-momentum shell of its arguments, and is inserted to take trace of the proper integration 
domains, that is 

M(fe}) = YlHm - A.) e (A - bll) . (106) 
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Up to this point, we assumed that the noise was zero mean Gaussian, white and TI. For sake of 
simphcity, we now assume that the noise represented by does not have spatial correlations as 
well, so that it satisfies (|73|). Eq.( |105| ) now reads 

d^q [q- {p- q)f M{q,p-q) 



-jd-p <t><i-p)\2D'X' 



>Mp)- (107) 



c) The third term of 0(A2) in Eq.([l02D adds to the original action a new term, not included 
previously (Fig. 3c.). The contribution of this term to the exponenciated CGA can be written as 



exp 



X^D / d\ 



An'^ (iq^ + uq' 

In turn, we can regard this contribution as coming from a new source of noise, in a sense that 
will become clear after we express e^^'^^th g^s the functional Fourier transform of an appropriate 



expression |3^, ^, ^ ^ . That is 

X exp\i{27Ty^ / d^qd^p p{q) (p<{-p)^<{p - q) 



{p- q) ■ q 



M(g) 



(109) 



{iq^ + z/g 2) 

where Z is normalization factor. Hence, e'^'^^"^ can be seen as the average of certain new term, 
according to the probability distribution of the auxiliary source p. This distribution is that of a 
white, TI and Gaussian noise, which has a second order momentum equal to 2DX^. Moreover, p is 
a multiplicative, rather than additive noise. 



d) The fourth term of order in Eq.( |102| ), is represented as a one loop diagram, built up by two 
propagators (0>v9>) (Fig 3d.). When calculating this propagator in the coordinate representation, 
((/)>(x)v9>(x')), the separation on lower and higher (spatial) wave numbers modes, does not prevent 
the arising of Si9{t' — t), as result of integrating p^ in the complex plane when the Fourier transform 
of {(p>{p)^p>{p')) is performed in reverse. Thus, the double product (0>(p')v'>(p)) (0>(p)v'>(p')) 
has null-measure support, and the diagram vanishes. 

(The fifth term is proportional to the propagator of two fields 0>, which is zero.) 



e) The sixth term is the sum of two non equivalent diagrams (Fig. 3e.): 

^2 



^— J d^pid^p2 0<i^<2 



+ 



j d'^psd^qid'^q2d'^q3 (27r) ^ 5123(5123 P2 " ^2 ■ % 
25(gi + g2) N{p3,q3) 



' {iql + uqi) {ipl + ypi){iql + uqi 

Niq2,q3) iSjqi+ps) 
{iq^ + uqi){iql + uqi) {ipl + upi 



M(p3,gi,g2,g3) 



(110) 
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The first delta function in tlie square brackets, gives a ^(gs), but % must be integrated in a sliell 
tfiat does not include the origin, and therefore the contribution of this member vanishes. For similar 
reasons, because N{p,p') oc Slp + p'), the second member in the square brackets also vanishes. 

f) The seventh term of order in Eq. ( |102| ) is represented by a single diagram (Fig. 3f.), which 
includes a loop formed by the propagator (0>iV>2) oc 6{qi + ^2)- This delta function and that of 
conservation 5j235 generate a S^qs). Because of the integration domain, as before, the diagram gives 
no contributions. 



g) The eighth term (Fig. 3g.) gives the following contribution to the CGA 

-A^ J d^Pi 0<(-pi) X (^{271)-^ J d%d% P3 ■ (pi-ps) gs ■ {Ps - Qs) (HI) 

Nipi -P3,P3 - qs) (p<iq3) 



{ip'i + upi){i[pi - P3T + - P3?)ii[P3 - 93]° + i^[P3 - gs] 



M(p3,pi -P3,P3 - ga) 



This contribution can be thought as a momentum dependent correction to the viscosity term in 
(|93|). When expanded in powers of the external momentum pi, the curly bracket takes the form of 
an infinite sum of derivative interactions. We saw in Section 5 that, if the noise satisfies (|7^), in the 
limit in which the shell is made of infinitesimal thickness, the expression between curly brackets 
gives, for small external pi, a factor proportional to p^. All other contributions are of higher order 
in \pi\. 

h) Finally, the ninth term of order A^ in Eq.( |102|) , generates a new vertex (the cubic interaction 
term), which couples three (p<s with one 0< (Fig. 3h.): 

- - J d p,d p,d q,d qs {2n) 5,,,, [-^ipU pl) + v{n + ^ ^<^^<^^<^^<^ ' (l^^) 

10 Appendix C 

In this appendix we compare the results of the previous section with those obtained by coarse 
graining the equations of motion. This is the first step of the transformation associated with the 
dynamical RG as defined in 0, to be further discussed below. 

10.1 Definitions 

In general, we start from a given stochastic equation 

c{^}{p) + Arw}{p) = v{p), (113) 

where the operator C is linear and Af collects the non-linear terms, and where, to be specific, 
the noise t] verifies ([7^). The non linearity couples modes of different scales. An exact solution 
can be attained in few cases only, such as the noiseless Burger's equation in 1+1 dimensions P8[] . 
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One could be interested in reducing the number of modes -for a computational calculation on a 



discrete lattice ^5] , or in study the scaling properties -in relation with critical phenomena 0. 
In both cases the elimination of short scale modes can be accomplished by solving their equations 
of motion in terms of the long scale modes, adopting, in general, some perturbative scheme. One 
then feedbacks these solutions in the equations of motion of the long scale modes, obtaining a 
coupled set of effective equations for these modes only. One identifies, in these equations, effective 
couplings -some of which were zero in the initial equations-, and noise terms, which can be either 
additive or multiplicative. 

Formally, we can define a projector V over the Fourier space spanned by modes in the momentum 
shell Kg <\p\ < A, and project the Eq.( |113| ) to obtain 



C{v>} + VM{'^y + V<} = V>, (114) 
£{(/.<} + (l-P)Ar{(^> + ^<} = v<- (115) 

In some way we must solve the first equation for ip^ to obtain y9>[(y9<, The second equation is 
then rewritten as 

+ (1 - V)Af{ip^ + ^>[ip<,v>]} = V<- (116) 

This will be the effective equation for the long modes, and the one we expect that reproduces the 
results obtained from coarse graining the OTP generating functional. Some fluctuating terms in the 
left hand side of ( |116| ) can be added to the noise rj^ to form an effective noise fj^, which will have 
an amplitude (or more precisely, a two point correlation function characterized by an amplitude) 
D. We remark that in Eq. ( |116| ) there is not implicit any kind of averaging process. The effective 
noise amplitude can be obtained trivially from ( 116 ) by calculating the correlation of f/<. 



This situation is different to that addressed, for example, in the paper of Medina et al. 
concerning the KPZ equation, where the effective noise amplitude is derived from the two point 
correlation function of the fields, or that presented by McComb for the Navier-Stokes equation in 
ref. EH], where the effective equation is averaged with respect to the short scale noise. For example. 



the right hand side of equation 9.34 in McComb's book p5| displays the unrenormalized external 



force, while our approach would replace it by the effective one (see eqs. 3.11 and 3.18 of 144|). This 



difference arises because of the way the average is performed in equation 9.16 of [45 



We show below the results of applying the coarse graining procedure to the KPZ equation. 
10.2 Coarse graining of KPZ equation 

Starting from Eq. (|93[) , we proceed as before, splitting the field as the sum of two independent fields, 
= V^> + 9?<, and analogously for the noise rj. Thus, 

{ipo + vp'^)[p> + </'<](p) + n (P- Q){f>i(l)f>iP - 9) + 2</?>(g)^<(p - q) (117) 



2 J 47r2 

+ <^<(g)<^<(p -q)} = [v>+ v<]{p)- 
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If IpI > As, then V5<(p) and ri^{p) are zero, and we obtain 

A f d^q 



{ipo + Jyp^)v>{p) + - J -^Q- (P- Q){v>{<l)v>{P -<l) + '^^>{(l)v<{P - q) (118) 

+ '^<{q)v<{p -q)} = v>{p)- 



This equation can be solved, formally, order by order in A by setting 

V.> = ^>W + Av^>« + ... 

It yields 

{ip^ + up'^) ' 

-1 r d^q 



v>''Kp) 



(119) 



(120) 



2{ip° + z/p2) J 47r2 

r]y{p-q) 



q- {p- q) 



V>{p-q) 



+ 2v;<(g) 



[i[p - g]° + i^[p- g]2) 



{iq^ + uq"^) {i[p — q]^ + ulp — g]^) 

+ ^<iq)v<ip - q) \ ■ 



(121) 



On the other hand, when \p\ < A^, fy{p) and f]>{p) are zero. For such p, and using the expressions 
given in eqs. ( 120 ) and ( |121| ), we find a closed equation for the field namely 



{ip° + iyp^)(p<{p) + ^ / ^ ^' (P~ ^ 

+ 2^<{p-q) 



^<{q)^<{p - q) 

V>{(l) 



+ 



v>{(i)v>{p - q) 



{iq^ + z/g 2) (^gO _|_ iyq'2'j(^i^p _ g]o _|_ iy[p — g]2) 



A^ 
2 



d^qd'^k 
167r4 



g ■ {p — q) k ■ (g — k) M(g) 



V<{p-(l) V<{k) <^<(g - fc) 



+ 2 



(zgO + z/g 2) 



+ 



{iqO + uq^){tk° + uk^] 
^<{p - q) V>{k) V>iQ - k) V <{k) ^ <{q - k) r]>{p - q) 



(122) 



+ 



{iq^ + vq^){ik^ + vk^){i[q-kf + v[q-kf) ' (ig° + z/g 2) (z[p - g]0 + zy[p - g]^ 

V9<(g - k) r]>{k) r/>(p- g) 



+ 2 



(ig*^ + vq'^){i[p — g]*^ + z/[p — g]2)(iA;0 -|- 

^>(^) v>{(i - k) v>{p-<i) 



{iqO + uq^){i[p - g]o + u[p- q]^){ikO + uk^){t[q - k]^ + u\q - kf) _ 
This is the basic result: an effective equation which only contains the long modes 



Tl<{p)- 



We now resort things in order to clarify the meaning of each term: The first term of 0{X) is 
the original non linearity. In the second one, ?7> acts like a multiplicative noise over and can 



29 



be identified with that we found earher in Eq.( |109| ). Rewrite the third term of 0{X) in Eg. ( |122| ) as 



A fd^q 



2 J 47r2 



q- {p- q) 



^ [ V>iQ)v>ip- q) - N{q,p- q): 

{iq^ + vq'^){i[p — q\^ + v[p — q\'^) \ {iq^ + i'q'^){i[p — g]° + z/[p — g]" 



Niq,p-q): 



.(123) 



With Ny we have indicated that the functions are zero if their arguments hes outside the mo- 
mentum shelL In this expression, the first term gives a field independent contribution, which when 
the noise is delta-correlated reproduces the result shown by Eq. (|104|) . In its turn, the remaining 
term in Eq.( p.23| ) is an additive source of noise, and therefore the effective noise term, to 0{X), is 
given by 



V>{(1)V>{P - g) - N{q,p- g)> 
+ h'q"^){i[p — g]° + ulp — g]^) 



(124) 



This noise has zero mean, and its two point correlation function, assuming i] satisfies (|7^, is given 

by 



2D + 2D'^\^ 



[q - {p-q)f M{q,p-q) 



A 

4^ [OT + Z/2(g'2)2][(p0_^0)2 + j,2{(p'_^2}2 



(125) 



Notice that the correction to the noise two point correlation introduced above is equal to that given 
in Eq.( p.07| ). However, the third order correlation function is not zero, so we cannot say that the 
effective noise is Gaussian, as was the original one. [However, the third order momentum is 0{X^)]. 
To say something about the higher order correlation functions first we must include corrections to 
the noise coming from higher order terms in the perturbative expansion we performed to arrive to 
Eq.O. 



Consider the O(A^) terms in Eq.(122). The first one is just the cubic interaction given in 
Eq.( |112|) . The second and the fourth terms are quadratic interactions subjected to multiplicative 
noise. The third term contains also multiplicative noise which, when the noise rj is TI, has zero 
mean (because of M(g), and, hence, it does not contribute to the effective viscosity. In the CGA 
all these multiplicative noisy terms appear when the perturbative expansion is extended to (9(A^). 



The fifth term (9(A^) can be written as 
A^ rd'^qd^k 



q - {p — q) k ■ {q — k) 



2^^{q-k) N{k,p-q): 



{iq^ + ijq'^){i[p — g]° + z/[p — q\^){ik^ + uk ' 



+ 



2v?<(g - k) {r]>{k) ri>{p - g) - N{k,p- g)>} 
{iq^ + vq'^){i[p — g]" + z/[p — q\'^){ik^ + vk ^) 



(126) 



Hence, we can regard the first contribution as the momentum dependent correction to the viscosity 
we found in Eq.( |lll[ ), and the second one as another term linear in subjected to multiplicative 
noise. This term, as the last one appearing in Eq.( |122[ ), which contributes to the effective noise, 
are found at O(A^) in the CGA. We conclude that the coarse grained equation of motion coincides 
with the equation of motion derived from the CGA. 
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